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A thermodynamical treatment of a massless scalar field (a "photon") confined to a fractal spatial 
manifold leads to an equation of state relating pressure to internal energy, PV S — U /d s , where d s 
is the spectral dimension and V s defines the "spectral volume". For regular manifolds, V 3 coincides 
with the usual geometric spatial volume, but on a fractal this is not necessarily the case. This is 
further evidence that on a fractal, momentum space can have a different dimension than position 
space. Our analysis also provides a natural definition of the vacuum (Casimir) energy of a fractal. 
We suggest ways that these unusual properties might be probed experimentally. 
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Among one of the more unusual properties of spatial 
fractals is that in some not yet completely understood 
sense the dimension of momentum space can be different 
from the dimension of position space. Even the notion 
of a (Fourier) mode decomposition on a fractal is not 
yet understood mathematically [D-Q- At first sight, this 
would appear to render problematical the conventional 
formulation of thermodynamics and statistical mechan- 
ics in terms of phase space cells Indeed, the quan- 
tization of the electromagnetic field and resulting notion 
of a photon appears inherently based on the Fourier de- 
composition of the field as a superposition of an infinite 
number of modes and their quantization as bosonic quan- 
tum harmonic oscillators. The very idea of mode decom- 
position is based on the assumption that it is possible 
to label and count the modes by quantizing momentum 
space into elementary cells of volume L~ d , where the spa- 
tial volume V — L d is taken to be a large hypercube in 
d space dimensions. 

While such a mode-based approach is not available at 
present for fractals, some results are known on fractals 
concerning the spectrum of the Laplacian, encoded in 
spectral functions such as the heat kernel trace and zeta 
function [l], 0, 0-0] . We argue this is sufficient informa- 
tion to define a consistent thermodynamic limit, lead- 
ing to an equation of state PV S = U/d s , in terms of a 
" spectral volume" V s [defined below] and the spectral di- 
mension, d s Q. On a regular manifold, the spectral and 
geometric volume coincide, as do d s and the geometric 
dimension d, but on a fractal this is not necessarily true. 

The idea of using thermodynamics to probe the ge- 
ometrical structure of objects goes back to H. Lorentz, 
who in 1910 addressed the original "Can one hear the 
shape of a drum?" question [9(, asking why the Jeans 
law, dE = V(T/TT 2 c 3 )uj 2 du!, only depends on the volume 
occupied by the black-body, and not on its shape. Soon 
after, H. Weyl found the geometric meaning of the large 
eigenvalues of differential operators , HH , and it turns 
out that this is what is probed in the thermodynamic 
large volume and high temperature limit. Physically, at 



thermodynamical equilibrium the photons probe the ge- 
ometric structure of the spatial manifold in which they 
are confined; the first geometric property we learn is the 
volume. We are also motivated by suggestions of fractal 
structure in models of quantum gravity 12 1. 

Our main technical result is an expression for the par- 
tition function Z(T, V s ) of quantum radiation at thermo- 
dynamic equilibrium at temperature T, on a fractal: 



In Z = y Jp Cr (In y) 



y = V s /(phc) d ° . (1) 



Here, (3 = l/(fesT), ks is Boltzmann's constant, and 
F pev is a periodic function of In (y), whose period depends 
on the characteristics of the fractal. 

We first review the usual relation between the parti- 
tion function and the effective Lagrangian for a mass- 
less scalar field in Euclidean space-time, emphasizing the 
importance of the asymptotic limit of large volumes or 
correspondingly lar ge t emperatures, and its relation to 
the Weyl expansion [13| . In this approach, the mode de- 
composition is encoded in the heat kernel trace for the 
Laplace operator, so we can then generalize to fractals. 

The partition function can be written as: 
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M is a circle of circumference h/3 in the coordinate t, A 
is (minus) the Laplacian defined on the spatial manifold 
A4, and c is the speed of light. While © is independent 
of any specific representation, it is usually derived using 
a mode decomposition of the field. Recall that the spec- 
tral partition function for a single bosonic oscillator of 
frequency lu is 
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By standard manipulations [22(, we obtain 
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Using the Kubo, Martin and Schwingcr (KMS) condition 
[Tij for thermodynamic equilibrium of a bosonic field <j) at 
temperature T, we recognize the second factor in the in- 
tegrand as a sum over Matsubara frequencies (jqlf) 2 , cor- 
responding to the spectrum of the operator d 2 /dt 2 with 
periodic boundary conditions <fi(t + hf3) = 4>(t). Idcnti- 
fiying ui 2 = c 2 k 2 with the eigenvalues of c 2 A, and tracing 
over all modes [23| , we recover expression ([2]) : 
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All information about the spatial manifold A4, including 
possible dependence on its volume V, is contained in the 
heat kernel trace factor. It is convenient to rescale r: 
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defining the "photon" thermal wavelength, Lp = (3hc, 
the dimensionlcss function f (t) = y^ + °° _ e~( 27m ' r , 
and ICm a s the dimcnsionless heat kernel trace of the 
Laplace operator on the manifold Ai 
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The heat kernel contains an implicit length scale L, char- 
acteristic of the geometry of the manifold Ai, that en- 
ables us to define a dimensionless Laplacian A = L 2 A 
[recall that the eigenvalues of A have dimensions of 
1/length 2 ]. Results in the mathematical literature refer 
to the dimensionlcss Laplacian A. 

As an illustration, consider black-body radiation in a 
large hypercube of volume V = L d , in d space dimen- 
sions. The corresponding heat kernel /Cy(r) is obtained 
from the mode decomposition of the field, u = c|k| = 
c27r|n|/V rl / rf , where n is a d-dimension integer-valued 
vector obtained by counting the elementary momentum- 
space cells of volume (2ir) d /V: 
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It is clear that the partition function has the form 

]nZ{T,V) = Q(x) , x = L p /V 1 ' d (9) 

where the precise form of the function Q is determined 
by the heat kernel trace in ([7|). Standard thermodynamic 
quantities and relations follow immediately from © . For 
example, the equation of state, PV = U/d, relating the 
internal energy U of the radiation to the pressure P and 
the volume V , follows immediately from 
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The Stcfan-Boltzmann expression for the internal energy 
U is a consequence of the equation of state and the ther- 
modynamic relation (f^) T = T — P, while notic- 
ing from (fTTj) that P depends on T only, in the thermo- 
dynamic limit. We then obtain U = aVT d+1 where a is 
a constant to be determined from ©. Also, the entropy, 

S = k B U/T + k B lnZ{T,V) =dxQ'(x)+Q(x) , (12) 

depends only on the dimensionless parameter x, so that 
an adiabatic (isentropic) expansion of the radiation oc- 
curs at constant values of the product VT d . 

For black-body radiation in manifolds of complicated 
shape, it is difficult to make an explicit mode decom- 
position and find an explicit expression like © for the 
heat kernel trace. However, ([9]) continues to hold. Wc 
can learn about the thermodynamic [large volume] limit 
from the Weyl expansion of the heat kernel trace, de- 
scribed below. Note that the large volume limit corre- 
sponds to V 3> Lp, which is a high temperature limit 
k B T hc/V 1 / d . This probes the small r behavior 
of 1Cm{ t ) for which a general asymptotic expansion is 
known: 
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where V is the radiation cavity volume, S is the surface 
area, and a is a numerical constant depending on the 
boundary conditions. Higher-order terms in the expan- 
sion characterize other geometric and topological prop- 
erties such as the curvature of the surface 
ing only the dominant volume term in 
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([6]) leads immediately to the familiar thermodynamic ex- 
pressions 0: ]nZ= {V/L d )C, R {d + 1)T (^±i) /vr^ 1 )/ 2 , 
P = (fc B T/L^)Cfl(d + l)r (^±i) l^ d+v >l 2 . Away from 
the thermodynamic limit, subdominant terms in f| 13[) 
lead to corrections that depend on the exact geometry 
of the volume enclosing the radiation, but the equation 
of state PV = U/d is always valid. This formulation in 
terms of heat kernels makes it clear that the Weyl expan- 
sion is directly related to the thermodynamic limit of a 
black-body radiation system, so wc can use the leading 
term as a definition of the volume probed by the photons 
as they attain thermal equilibrium. 

We now turn to fractals, and make use of recent 

S ogress concerning the heat kernel trace on fractals 
Q, H-Q- While it is familiar that fractals have a 
nontrivial fractal (Hausdorff) dimension dh, arising from 
their self-similar spatial structure, it is perhaps less well- 
known that fractals are characterized by another dimen- 
sion, known as the "spectral dimension" d s , which need 
not be equal to dh [8(. Indeed, diffusion on fractals is 
"anomalous" in the sense that the Einstein Brownian 
motion relation, (r 2 (t)) oc t, for long enough times, is re- 
placed on a fractal by (r 2 (t)) oc t 2 / dw , where d w = 2dh / d s 
is called the "anomalous random walk dimension" |8f. 



3 



K(t) ^eadi„g (t) 



0.8 
0.6 
0.4 
0.2 




0.00 0.05 0.10 0.15 0.20 0.25 



by the appearance of oscillatory terms arising from the 
behavior of the log periodic function F in (|14l) . 

In order to make this more explicit, and to give the ac- 
tual numerical coefficient appearing in the spectral vol- 
ume, it is convenient to make a Mellin transform to con- 
vert our expressions from the heat kernel trace (|7|) to the 
associated (dimensionless) zeta function : 
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with inverse (C is the usual inverse Mellin contour) 



FIG. 1: Typical small r behavior of a fractal heat kernel trace. 
There are small oscillations [see inset], as a function of lnr, 
about the leading behavior ~ l/r da ' /2 . 



For many fractals, including well-studied cases such as 
the Sicrpinski gaskets or diamonds, the small r limit 
of the heat kernel trace is more complicated than the 
Weyl expansion expression in (|13[) . The heat kernel trace 
K.jt{t), on a fractal J 7 , for the dimensionless Laplacian 
A, behaves as 1Cf(t) = T~ d ' B ^ 2 F(t), where F is a pe- 
riodic function of lnr The typical small r behav- 
ior is shown in Figure 1. The leading small r behavior, 
^■J r i T ) ~ T - d B/2^ j g determined by the spectral dimen- 
sion d s , but there are also oscillations around this lead- 
ing behavior, and these oscillations are typically very 
small. Explicit expressions for this oscillatory behav- 
ior are given in Q for diamond fractals, and numerical 
results for the Sierpinski gasket appear in (3). Restor- 
ing units for the Laplacian, there exists a characteristic 
length L s (a "spectral length") such that, 
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Inserting this expression into ^ leads to (U} for the par- 
tition function In Z. 

The first important conclusion from this result is that 
the thermodynamic volume is determined by the spectral 
dimension d s rather than the fractal dimension dh- So, 
we define the thermodynamic volume of the fractal to be 
the "spectral volume" : V s = b L ds , where b is a numerical 
coefficient that we specify in detail below. This implies 
that the thermodynamic equation of state on a fractal is 



PV a = 

d s 
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and according to (|12p , adiabatic (isentropic) transforma- 
tion on fractals occur at constant values of the product 
V s T ds . The second important conclusion is that the ac- 
tual expressions for pressure P, internal energy U, etc... 
will be modified on a fractal, not only by the appearance 
of the spectral dimension and spectral volume, but also 



K, m {L}t) = — 
H 2m 



ds 



c 



T^) Qm{s)T{ S )t- s . (17) 



Straightforward manipulations show that the partition 
function can now be expressed as 
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The first term gives the standard zero temperature "vac- 
uum energy" contribution [l3j], proportional to £m (— f)j 
which has recently been generalized to quantum graphs 
17j . On a fractal J 7 , our thermodynamical analysis gives 
a simple expression for the Casimir energy: 
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The second term in (|18[) encodes finite temperature cor- 
rections to the internal energy U ; the Ricmann zeta fac- 
tor £fj(2s+l) arises from the sum over Matsubara modes. 

On a regular manifold M. , the zeta function C,m (s) is a 
mcromorphic function in the complex plane with simple 
poles on the real axis. The pole with largest real part 
is at s = d/2, and gives the leading term in the Weyl 
expansion (|13j) . Thus, we find a natural thermodynam- 
ical definition of the volume in terms of the residue of 
Cm(s) at this pole: V = (4ir) d / 2 T(d/2)L d Rcs(d/2). De- 
pending on the form of the manifold, and the boundary 
conditions, there may be other poles on the real axis, and 
these determine the sublcading terms in (TIB")) . When the 
spatial manifold M. is a fractal, the situation is radically 
different 0, 0, [3, [lj| : the zeta function of the (dimen- 
sionless) Laplacian A has also been shown to be mero- 
morphic in the complex s plane, with simple poles, but 
now the poles can be complex, as illustrated in Figure 2. 
The pole with the largest real part lies on the real axis 
at s = dg/2, where d s is the spectral dimension. This 
pole leads to the dominant behavior, shown in the pre- 
factor in (| L4[) , and provides the definition of the spectral 
volume: 
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FIG. 2: The generic pole structure of the zeta function on 
a symmetric fractal [la . [T^ | . The dominant behavior comes 
from the real pole at s = d s /2, and oscillations from the tower 
of complex poles in (|21|) . There may be other poles in the 
shaded region to the left of this tower, but these do not con- 
tribute to the thermodynamic limit. On the Sierpinksi gasket, 
there is another tower of complex poles on the imaginary axis. 



In addition, there is an infinite tower of complex poles 
with the same real part: 
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is split into in each iteration [6, H, 19j|. These complex 
poles are the origin of the log periodic oscillations in (fT4"|) , 
and have been identified with complex dimensions for 
fractals [H, [H| . Explicit expressions for the zeta function 
on diamond fractals lead to simple expressions for these 
poles and their residues Q . Integral representations have 
been given in (l9j for the zeta function on other fractals, 
including the Sierpinski gaskets. In Q it was argued 
that the origin of these complex poles is the appearance 
of exponential degeneracy factors, rather than power-law 
degeneracy factors that appear for regular manifolds. 

To conclude, we have shown that the thermodynamical 
volume of a fractal is determined by the spectral dimen- 
sion and the "spectral volume" defined in (f20)). and fur- 
thermore that thermodynamical relations will have small 
oscillatory behavior in the thermodynamic limit. These 
two properties serve as direct and clear evidence of an 
underlying fractal structure. To observe such phenom- 
ena, we propose studies of thcrmalization of photons in 
fractal-shaped (e.g., Sierpinski) mesoscopic waveguides. 
Of course, one cannot fabricate such a waveguide to all 
orders of fractal iteration, and the photon wavelength 
sets a natural cutoff scale. However, computations show 
that the role of the spectral dimension, and even the os- 
cillatory behavior from the complex poles, can be seen af- 
ter only about five iterations, which should be accessible. 
A similar argument appears in [2(3]. Another interest- 
ing example is provided by systems such as a Fibonacci 
distribution of dielectric layers [2l[ , whose spatial struc- 



ture is not fractal, but whose energy spectrum is given 
by a fractal triadic Cantor set. The study of such lay- 
ered structures would be significant, because here the ge- 
ometric volume is well defined, and non-fractal; but we 
predict that the thermodynamically significant volume is 
the spectral volume. 
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